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Abstract 

If light supersymmetric top (stop) quarks are produced at the LHC and decay via on- or off-shell IK-bosons they 
can be expected to contribute to a precision W~^W~ cross section measurement. Using the latest results of the CMS 
experiment, we revisit constraints on the stop quark production and find that this measurement can exclude portions of 
the parameter space not probed by dedicated searches. In particular we can exclude light top squarks up to 230 GeV 
along the line separating three- and four-body decays, h —)■ x\W^*''b. We also study the exclusion limits in case when 
the branching ratio for these decays is reduced and show significant improvement over previously existing limits. 
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1. Introduction 

Searches for stops — the supersymmetric (SUSY) part¬ 
ners of top quarks — have received significant attention 
from both ATLAS [1-6] and CMS [7-12]. While limits 
obtained after Run 1 of the LHC at -^/s = 8 TeV can go, 
depending on the decay modes studied, up to 800 GeV, 
there are still parts of parameter space where relatively 
light stops are allowed, see e.g. the summary plots by AT¬ 
LAS [13] and CMS [14]. 

The main motivation for light stops is the so-called nat¬ 
ural supersymmetry [15] paradigm which demands that 
the particles must be close in mass to the ordinary top 
quark. Unfortunately however, this region of parame¬ 
ter space is particularly difficult to explore due to the 
background of top quark production. In particular if the 
stop quark decays via a top quark that is almost on-shell 
(U ^ no exclusion limit is currently present. An¬ 

other difficult region of the parameters space can be iden¬ 
tified at the border between three- and four-body decays 
with a (nearly) on-shell W boson (U —)• XiW^*^b). 

Several recent theoretical studies have attempted to 
fill these holes in the stop parameter space by using pre¬ 
cise predictions and measurements of top quark cross sec¬ 
tion [16] (see however ref. [17] for a discussion of pos¬ 
sible problems with this approach), specialized mono-jet 
searches [18], recasting other SUSY searches [19] or via an¬ 
gular correlations [20]. A complementary idea is that cer¬ 
tain corners of the parameter space might be constrained 
by looking for signals of stoponium production [21]. 
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An alternative approach presented here is based on 
the observation that light stops decaying into certain final 
states can contribute to the W~^W~ cross section mea¬ 
surements [22-26]. Until recently the ATLAS and CMS re¬ 
sults were displaying a moderate excess over the standard 
model (SM) prediction [27-29] but this was determined 
to be the result of neglected higher order corrections [30- 
32]. In any case, the fact that the observed cross-section 
was greater than the predicted background meant that 
any derived constraints on stop production would have 
been weak. However, the recent CMS measurement [33] 
based on the full Cs = 8 TeV dataset, using the next- 
to-next-to-leading-order (NNLO) cross section prediction, 
_)• 1U+1U“) = 59.8^];^ pb [30], and event re¬ 
weighing [32] turned ont to be very well aligned with the 
SM: cr®’'P = 60.1 ± 4.8 pb. In this Letter, we recast the 
CMS analysis as a potential way to constrain the produc¬ 
tion of light stops. 

We focus on three widely studied decay modes that are 
commonly present in SUSY models with light stops and 
improve the existing constraints. Assuming that only the 
light stop and the lightest supersymmetric particle (LSP, 
in our case the lightest neutralino, x)*) have masses of order 
of the electroweak symmetry breaking (EWSB) scale we 
have: 




if 

> wt -1- m^o , 

(1) 


^xlwb, 

if 

Hiii ^ ''^w + mb + m^o , 

(2) 



if 

< mw + mb + m^o . 

( 3 ) 


The three- and four-body decays might compete with loop- 
mediated two-body decay, ii —)■ x^c [34, 35], but the 
branching ratios (BR) are highly model dependent here [36- 
39]. Another possibility is given by: 

ii^Xib^XiW^*H, (4) 
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provided m-± < mz . Depending on the parameter point 
under consideration, in particular the mass differences be¬ 
tween stop and electroweakinos and their mixing charac¬ 
ter, the chargino mediated (4) and one of the direct de¬ 
cays (1), (2) or (3), may be simultaneously present, see 
e.g. [40]. This feature would have a significant impact on 
the expected exclusion limits. 

2. Simulation 

Monte-Carlo stop samples with up to one additional jet 
were simulated using MadGraph5_aMC@NL0 [41] and matched 
to the Pythia 6 [42] parton shower. The cross sections 
were normalized to the next-to-leading-order (NLO) pre¬ 
diction using NLLfast [43, 44]. 

The analysis of the simulated samples was performed 
using CheckMATE [45, 46] and a dedicated implementation 
of the CMS W+W- cross-section measurement. CheckMATE 
uses a specially tuned version of the Delphes 3 detector 
simulation [47] and jets were clustered using Fast Jet [48] 
with the anti-Zc-r algorithm [49]. The analysis is performed 
for di-lepton final states with missing energy which for the 
signal process originates from neutrinos. In order to sup¬ 
press a dominant SM background, ti production, events 
with 6-jets and multiple final-state jets are vetoed. To bet¬ 
ter understand the tt background, CMS defines two signal 
regions (SR): 0-jet SR without jets with pT > 30 GeV; 1- 
jet SR with exactly one jet with pr > 30 GeV. These are 
further subdivided based on whether the final state lep¬ 
tons have different (e/i) or same flavour (ee or pp). The 
expected and observed event numbers agree well within 
errors for all SR and therefore the analysis can serve as a 
constraint for models that contribute to the similar final 
state. 

Our implementation was validated using the event num¬ 
bers provided by the CMS collaboration for W^W~ sig¬ 
nal and tt background [33]. The samples used for vali¬ 
dation were obtained using MadCraphB.aiMCONLO [41] and 
hadronised using Herwig++ 2.7 [50, 51]. The parton dis¬ 
tribution function (PDF) sets used were CTEQ6L [52] for 
leading order generation and CTIO [53] for NLO. 

In order to apply limits to stop production we imple¬ 
ment two different procedures. For the first we calculate 
the model independent confidence limits in a modified fre- 
quentist approach (CLs method [54]) at 95% for the two 
flavour inclusive SRs. We then define the stop model as 
excluded when it predicts a cross section in excess of any of 
these limits. The second method produces a more strin¬ 
gent model dependent CLs limit at 95% by performing 
a combined fit to the four separate signal regions. We 
use HistFitter [55] as an interface to HistFactory [56], 
RooStats [57], RooFit [58] and ROOT [59, 60], to allow all 
background sources to float independently whilst including 
the correct correlated systematics. 

Many of the backgrounds to W^W~ production are 
determined with a data driven technique in the CMS anal¬ 


ysis. Essentially control regions are defined for the vari¬ 
ous backgrounds that are only expected to contain a small 
contribution from the W~^W~ process under study. These 
are used to normalise the various backgrounds and Monte- 
Carlo is then used to extrapolate to the W~^W~ signal 
regions. One may therefore worry that stop production 
contributes in these control regions, spoiling the normali¬ 
sation constants. 

Unfortunately CMS does not publish the control re¬ 
gions used so it is impossible for us to explore these effects. 
However, we note that it is possible that stop production 
contributes to the control region and thus increases the 
value of the normalisation constant. In turn this could 
increase the background prediction in the W~^W~ signal 
regions. Such an effect would actually strengthen the limit 
we derive on stop production and one may worry about set¬ 
ting a spurious exclusion but we believe that such a pos¬ 
sibility does not exist in our study. Firstly, the predicted 
and measured W~^W~ cross section are now in very good 
agreement suggesting that any signal contamination can 
only be slight. Secondly, of most concern for this analysis 
is the tt control region. However, in the parameter region 
of most interest for our study, nif « mi, + mw + , this 

corresponds to the 6-quarks being extremely soft. Hence 
it is very unlikely, that this final state will contribute to a 
tt measurement at all. 

3. Results 

The obtained exclusion limits for stop decaying via (1), 
(2) or (3) with 100% branching ratio are shown in Fig. 1. 
The best exclusion limit is obtained along the line ~ 
mw+iTib+m^o and is exactly where no current LHC search 
sets a limit on these models. The W~^W~ measurement 
allows us to constrain stop masses up to ~ 220 GeV for a 
LSP mass of ^ 130 GeV. We also note the additional exclu¬ 
sion for stop masses ~ 170-190 GeV for decay mode (2) 
where an intermediate top quark is nearly on-shell that 
also extends the limits from dedicated stop searches. 

The reason that the W~^W~ cross section measurement 
is so sensitive along this line is that the final states are most 
similar to the actual SM production of 1U+1U“ pairs: the 
W boson is (nearly) on-shell and the 6-jet is rather soft, 
significantly reducing 6-jet veto effectiveness. On the other 
hand, this region is problematic for dedicated stop searches 
due to its similarity to the SM background. Dedicated stop 
searches attempt to place cuts that act as a discriminator 
between signal and background. However in regions where 
the signal has very similar features to the background this 
approach breaks down and consequently our approach is 
complementary to other searches. 

It was shown in Ref. [39] that the branching ratio for 
decays (2) and (3) can be substantially reduced in favour 
of loop-mediated flavour-changing two-body decay, ti —> 
Xl c. Such a reduction can pose a significant challenge 
for dedicated stop searches as can be seen for example in 


2 




Figure 1: The exclusion limits for stop pair-production in plane assuming that only decay modes (l)-(3) are allowed, respectively. 

The dotted-blue line denotes the exclusion using 0-jet signal region and the red-dashed 1-jet signal region of Ref. [33]. The black solid-line is 
for the combined exclusion, as discussed in Sec. 2. The experimental exclusions were extracted from the following studies: ATLAS monojet [5], 
ATLAS 1-lepton [2], ATLAS 2-lepton [4], CMS 1-lepton [9], ATLAS spin correlation [6]. 


Fig. 12 of Ref. [9]. Therefore, we compare exclusion limits 
obtained in the current study with the results reported by 
collaborations, but now assume BR{ii —>■ = 0.6. 

In Fig. 2 we see that the limits are severely weak¬ 
ened and much more of the stop parameter space is un¬ 
constrained.^ The W^W~ measurement is still effective 
however and allows us to constrain stop masses up to 
~ 180 GeV for a LSP mass of ~ 80 GeV. In fact for 
low stop masses we can successfully exclude models with 
m^o < 60 GeV which current searches are not sensitive to. 

4. Summary 

We analysed constraints on the stop sector in light of 
the recent measurement of W~^W~ production cross sec¬ 
tion by GMS. We show that this measurement provides 
constraints on light top squarks that are complementary 
to the dedicated LHC searches. 

The best sensitivity is obtained along the line where 
an intermediate decay-mediating W boson becomes on- 
shell, where the conventional stop searches have a particu¬ 
lar weakness. Assuming that this the only available decay 
mode, the reach is mj:_^ > 230 GeV. The method retains 
its sensitivity even for significantly reduced branching ra¬ 
tios to the analysed final state. We demonstrate that in 


^The current available data for the ATLAS monojet study [5] do 
not allow for a reliable combination when branching ratios of less 
than 100% exist in a model. In addition, since the study relies on 
charm tagging that is difficult to reliably simulate with a fast detector 
simulation we remove the study from this figure. 


case of BR = 0.6, stops with masses < 180 GeV are 
also excluded. We note that in the reduced branching frac¬ 
tion scenario, the other searches are significantly limited 
and this additionally shows the complementarity of these 
approaches. 

Note added 

After completing this study, a summary of ATLAS 
Run-1 stop searches has been published [61]. It includes 
a dedicated stop search along similar lines to the sugges¬ 
tion in this Letter. We note that our results for 0-jet SR 
are consistent with those presented in ref. [61]. However, 
in some parts of the parameter space, in particular for 
~ TO(, m^o ~ 0, the GMS 1-jet SR [33] offers a stronger 
bound. 
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